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Since the CP violation in the B system has been investigated up to now only through processes 
related to the B-B mixing, urgently required is new way of study for the CP violation and estab- 
lishing its origin in the B system independent of the mixing process. In this work, we explore the 
exclusive B° — * p°7 decay to obtain the time-dependent CP asymmetry in b — *• d decay process in 
the standard model and the supersymmetric model. We find that the complex RL and RR mass 
insertion to the squark sector in the MSSM can lead to a large CP asymmetry in b — > d*y decay 
through the gluino-squark diagrams, which is not predicted in the Standard Model induced by the 
B-B mixing. 



I. INTRODUCTION 



The unitarity of the Cabibbo-Kobayashi-Maskawa (CKM) matrix yields the relation V ud V* b + V cd V* b + V td V t * b = 0, 
which describes a triangle on the complex plane. The angles of this triangle are defined by 

^ 2 )=arg(-##V /3( W = argf-§#V 7 (fc) - arg (-^) , (1) 
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of which non-zero values indicate the existence of the CP violation. In the B meson system, the CP asymmetry in 
B — > J/tpKs decays, which is directly related to sin 2/3 through B°-B° mixing, has been measured and agrees well 
with the standard model (SM) prediction by the CKM unitarity. Since the present CP violating asymmetry appears 
via /3, presently measured only through processes involving the B°-B° mixing, it is strongly required to find new 
observables of the CP violation in the B system and establish that the measured /3 indeed originates from Vt d of the 
CKM in a way independent of the mixing. Within the SM the source of the CP violation of the B system (and the 
weak phase /3) is nothing but the phase of CKM matrix element Vtd in the leading order. The magnitude and the 
phase of Vt d cannot be measured in a direct manner due to top quark's fast decay before its hadronization. This 
leaves only indirect methods for experimental determination in various processes. 

As is well known, sin 2/3 can be determined independently of the B°-B° mixing in a theoretically clean manner, up 
to A sin 2/3 w ±0.10, provided B(K + — > tt + vv) and BiK^ — > Ti a vv) are measured within ±10% accuracy 1]. Why do 
we urgently require yet another method to examine the CP violation in the B system? The main reason is that this 
area is very likely to yield information about new physics beyond the SM. We expect that new physics can influence 
the AB = 1 penguin decays in a different way than the usual AB = 2 mixing. Examples are the possibly large new 
physics effects recently measured in b — > sss decays. The measurements of sin 2(3 in b — > ss decays may have shown 
deviations from that in B — > J/ipKs decay, which are still controversial 0. A sizable difference between them could 
be a clear evidence of the new physics beyond the SM. If one of such discrepancy indeed exists, it implies an evidence 
of a new physics effect beyond the SM [3|, and also suggests an interesting possibility that the new physics effect is 
large in the penguin diagrams while its contribution to the B°-B° mixing is small. 

The Cabibbo-suppressed b — > d'j decay also involves Vt d and provides us a new chance to measure the weak phase 
/3 in a way independent from the mixing to study the CP violation in the B system. Although the inclusive B — > X d "f 
decay is a theoretically clean process to determine Vt d Q, it can be hardly discriminated from large B — > X s j 
background in the experiment. Thus our interest is devoted to the study of the exclusive channels, B —> pj and 
B — ► CJ7. The first measurement on the branching ratio of b — > process, 

Bt(B -> p/«7) = (1.8±°; 6 5 ± 0.1) x lO" 6 , 

has been reported at the Moriond conference [a. Implications of B — > p7 and B — * u^f decays in the SM and the 
MSSM have been studied in the literature H,0,|8(. The charged B ± —> p ± 7 decay mode provides clean signal and has 
a branching ratio twice larger than that of the neutral mode, by the isospin symmetry. However, the long-distance 
(LD) effect on the charged mode due to dominantly W ± -annihilation is very large (~30 %), which contaminates the 
CP violating effect [UEI3- Consequently here we study the B° — > p°7 decay to explore the CP violation. The time- 
dependent CP asymmetry in the neutral B° — > p°7 decay probes the CP violation in the interference between decay 
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and mixing as well as the direct CP violation, and it leads to the extraction of (5 from both decay and mixing. We 
note that the photon has two helicity states 7^ and 7^ which cannot be discriminated in the B factory experiments. 
Since CP violating asymmetry is defined when both B and B mesons decay into a same state, what we actually 
measure is the CP asymmetries with the definite helicity in the final states. In the SM, the operator which governs 
b — > dj decay is chiral and the term with the conjugate operator is suppressed by m d /m b and the CP asymmetry also 
suppressed accordingly. Therefore, the new physics beyond the SM which involves a sizable ri ght -handed operator is 
required for a large time-dependent CP asymmetry enough being observed in the experiment 

In this letter, we consider the supersymmetric models which have non-diagonal elements of the squark mass matrices 
in the so-called super-CKM basis as a new physics model. The non-diagonality of the down-type squark mass matrix 
is parameterized by the mass insertions 

(Sij)MN = {m%)MN/m 2 , 

where m is the averaged squark mass. Here i and j are flavor indices and M and N denote chiralities, L or R. 
The 5's are complex in general and provide new CP phases. To clarify and simplify our discussion, we consider the 
(Si3)rl and (Si3)rr dominating cases to produce the right-handed operators in the present work. We note that the 
($13)lr and (6i3)ll mass insertions leads to the left-handed operators and cannot generate the CP asymmetry only 
by themselves. The AB = 1 decay amplitude depends on (5i3)rl : rr linearly, while the dependence of M\ 2 for the 
AB = 2 process is quadratic. Then b — ► d'y decay would be more sensitive to the SUSY contribution than the B — B 
mixing process when the SUSY effects are smaller than the SM values as usual. Thus we can observe the new physics 
effect through decay process which is not observed through mixing induced processes. In section II, we describe the 
B° — > p°7 decay and the time-dependent CP asymmetry. The supersymmetric contributions are given in section III 
and the numerical results given in the section IV. We conclude in section V. 

II. CP ASYMMETRY IN 5° -> p° 7 DECAY 

The relevant terms of the effective Hamiltonian for the b — > g?7 decay is written as 

A(~* 

H *s = E [V&V^dOl + CzOl + C^ + c 2 o' 2 q ) , 

* q—u^c 

~v tb v; d (cfo 7 + c'fo' 7 ) + •••], (2) 

of = fea/xsfyio, 

0% = (d LW )(q L rb L ), 

O- = {em b /16ir 2 )d L <j^F^b R , (3) 

and 0[ are their complex conjugate operators. The effective Wilson coefficient c\'^ cS includes the effects of operator 
mixing. 

In the SM, C 7 eS is suppressed by the mass ratio m^fmh and so is the right polarized photon emission b^ — > q R jR. 
The photon emission with the momentum q and the polarization vector e M (g) yields the electromagnetic field strength 
Ffiv — i{tvqii — e /-iQv)- The hadronic matrix element of B — > pj process is parameterized in terms of two invariant 
form factors F£ and F| and the polarization vector of the p meson, e M (fc). We can set F[ — F 2 S = F s at q 2 ~ for 
the real photon emission ^3] ■ We write the amplitudes for the final states of polarized photon as 

A L = { PlL \H cS \B ) = ^C 7 cS *V t lV td (piL\O' 7 t \B"), 
AC 1 

A R = ( P jr\H cS \B q ) = -j£- C f*V t tV td (pjR.\OS\B°) 7 
A L ( PlL \H cS \B ) = ^CfV tb V; d (piL\0 7 \B°), 

A R EE { PlR \H cS \B°) = ^-C' 7 cS V tb V: d (p lR \O' 7 \B ), (4) 



where the operators are given by 
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where 



{p{k) lL (q)\0 7 \B\p)) = ^e»(qy(k)[e fl „ al3 p<*qe-i(g fl „(q-p)- Pfl q„)]-2F s (q 2 ~0), 

= (p lL \0'S\B°), 

(p(k) lR (q)\O> 7 \B (p)) = ^e»(qy(k) [e^ a pp a q^+i(g^(q-p)-p^)} ■ 2F s (q 2 ~ 0), 

= (pUOr^B"). (5) 

For the neutral B meson decay, the LD contribution due to W-exchange is merely a few % from the QCD sum rule 
calculation so it will be ignored in our analysis. We investigate the time-dependent CP asymmetry of neutral 

B mesons, which probes the direct CP violation and extracts the weak phase in mixing and decay. Since we cannot 
measure the photon helicity in practice, the asymmetry which we will measure in the experiment is given by 

, (f) = [r(^° hvs ^ p°il) + r(g° hvs -> - [r(i?° hys -> Al) + r(^ hys -> P ° lR )] 

CH ' nB° phys - P°1l) + T(B° hys - pO lR ) + T(B° hys - p° 7i ) + r(B° hys - p° 7fl ) 

= — Ccos(Amst) + <Ssin(Amst), (6) 



where 



\A L \i + \A R 

with the parameter A^(^) defined by 



_ \A L \ 2 Im\ L + \A R \ 2 Im\ R 

* n 12 ' UJ 



_ /M^A L(R) 

V Ml2 Ac(.R) 

Since |j4l| = and |Ar| = the coefficient C = identically. The off-diagonal element Mi 2 describes the 

B°-B° mixing and ^l(_r) is the amplitude for the b — > d 7 L(ii) decays. We define 

/3 mix = -Arg(Mi 2 ), 
/3dcc ay = ^Arg ( *L) = iArg f ^ , (9) 



2 2 °V^l 

and the coefficient S is expressed by 



5 : -^T^ sin(2/3mix " 2/3docay) - !lm 



We can write 

2/3dcca y = 2/3 SM + Arg(C£) + Arg(C 7 ). (11) 

Note that we have an additional factor 2 \Cj\ | C"^ | / ( | C7 1 2 + | | 2 ) , which can enhance or suppress iS by the new physics 
effect \C 7 \. In the SM, | Cy | / 1 Cr | ~ 0(m<2/m&) and the CP asymmetry is also suppressed by the corresponding factor. 

III. SUSY CONTRIBUTIONS 

We consider the gluino mediated penguin diagram contribution to b — > g? 7 decay in the MSSM. By penguin diagrams 
with gluino-squark loop, the Wilson coefficients C[ in the effective Hamiltonian of Eq. (2) get contribution to produce 
j R at the matching scale p = mw- After the RG evolution, we have 

Cf{m b ) = -0.31, 

C'f (m b ) = _ ^- (0.67 C^ SY (m w ) + 0.09 Cf USY (m w )) , (12) 
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where the SUSY contributions are 

4a s 7rQb 



Cl VSY (m w ) 



3rh 2 



(Si 3 )rrM 4 (x) - (6ia)aL4Bi(jc)- 



m b 



6rh 2 



(Sx 3 )rr{9M 3 {x) - M A {x)) + {6 13 )rl 4£?i(x) - 9 



B 2 (x)\ m 



(13) 



with the squark mass insertions (£13) bt, and (6i 3 )rr and the squared mass ratio of the gluino mass to the average 
squark mass x — (m g /fh) 2 [l3l [T3. Ilq. Note that the SUSY contribution is more sensitive ^ to ql than (Sx 3 )rr 



due to the enhancement factor rrig/mb- The loop functions Bi(x) are found in the literature |13lll4Lll5l |. Since Srl,rr 
are complex in general, the Wilson coefficients C' e 7 (nib) has nontrivial phase which affects the phase of A/ A. 

On the other hand, the B-B mixing is affected by the gluino-squark box diagrams in the MSSM. Unlike the AB = 1 
case, the relevant AB = 2 effective Hamiltonian with the supersymmetric contribution contains new operators which 
consist of scalar-scalar interactions, 



0' S 2 = (d*(l+76)U0*/s(l + 75)6/3), 
O'ss = (4(1+75)M(4(1+75)M, 



(14) 



when we introduce only the RL and RR mass insertions. The Wilson coefficients C\ t s2,S3 corresponding to the 
SM operators 0\ = [dy^l — ^)b)(d'y ll (l — 75)6) and Os2,s 3 — 0' S2 S3 {L <-» R) consist of the SM part and the 
supersymmetric contributions, while C' S2 and C' S3 corresponding to the above operators are entirely supersymmetric. 
Their explicit expression at the scale /i = -Msusy can be found in Refs. [Til Il7j . Ignoring the RG running effects 
between -Msusy and mw, we perform the RG evolution from mw to nib scale to obtain the evolved Wilson coefficients. 
The SUSY contributions are given by 



C K >i (m b ) 



^cW- usy (m ; 



SUSY 



), 



(15) 



(id) 



where the magic numbers a,k, b 
element M 12 obtained by M 12 = {B°\H A 
in vacuum insertion approximation. 



and c 
AS 

eff 



are found in Ref. 0] and i] — oi s (MsuSy)/ 'a s (mw)- The off-diagonal 
2 \B°) /2niB consists of the bag parameters Bi and the decay constant fs d 



IV. NUMERICAL RESULTS 



We write {Si 3 )rl(rr) = \($13) RL(RR)\e lip ■ Figure 1 shows the allowed values of the quantity S as a function of the 
phase ip, assuming (5i 3 )rl dominating case with |(<5i3)_rl| = 0.001. We vary the weak phase 7 from to 2n. Hereafter 
we use the input parameters as follows: nis — 5.3 GeV, m t — 174.3 GeV, m b — 4.6 GeV, and a s (mz) = 0.118. The 
decay constant /s d ~ 200 ±30 MeV is the main source of the theoretical uncertainty and the bag parameters are those 
of Ref. 0| B\ = 0.87, B 2 — 0.82, B 3 — 1.02. The supersymmetric scale is taken to be nig ~ rh ~ Msusy ~ 500 
GeV. We require that the mass difference Awib and /3 m ix in B —* J/jbK decay should be within the experimental 
limit: Am B = 0.489 ± 0.008 ps _1 ^ and sin2/3 mix = 0.734 ± 0.055 0. On the other hand, it is reported that the 
branching ratio of the inclusive B — > Xd"f decay puts a stringent constraint on the flavor mixing parameters of the 
gluino-squark contributions (l7|. In spite of the recent measurement of the decay rate of the exclusive B — > p/uj'f 
channel, we do not use Br(B — > p/tu^/) as a constraint here since it involves a large theoretical uncertainty in the 
form factor. Instead, we assume a moderate upper bound Br(B — * Xdj) < 1.0 x 10~ 5 in this analysis, following Ref. 
m and referring to the measurement of the exclusive channel 0, although no experimental limff on the branching 
ratio for the inclusive decay is given yet. In the figure, the black region corresponds to the allowed values for the 
phase of (S± 3 )rl, while the grey (green) region to the parameter sets which satisfy the Am# and sin2/3 m j x constraints 
but exceeds the bound on the branching ratio of the inclusive B — > Xd"f decay. We find that large CP violating 
asymmetry reaching ~ ±60%, is possible while all other experimental constraints are satisfied. As we can notice from 
Eq. (10), the maximal value of the CP violating asymmetry is set by the coefficient 2 |CV| | C*^ | / ( | C7 1 2 + |Cy| 2 ). In the 
present case, this coefficient is roughly 0.6 because the allowed value of \C' 7 \ is about 0.1 for |(5i3)_rl| = 0.001 (note 
Cj = Cf M = —0.31 for {Si 3 )rl(rr) dominating cases). If we increase the assumed value of |(<5i3)kl|, the allowed 
range for the CP asymmetry, S becomes expanded and may reach even ±100% if |CV| = \C' 7 \, which is realized 
when |(<5i3)_rl| ~ 0.003. However, such large \C' 7 \ values would give large B — > Xdj branching ratio which exceeds the 
assumed upper limit of the branching ratio. Still the CP asymmetry S ~ ±0.8 is allowed while satisfying experimental 
constraints including the inclusive branching ratio bound. 
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The plot for the allowed range of S with respect to |(<$i3).rl| is depicted in Fig. 2 when the phase ip is fixed to be 
zero. The black region and the grey (green) region are defined as in Fig. 1. We see that the value of the CP asymmetry 
is zero for |(<5i3).ri,| = 0, which corresponds to the SM case with vanishing d- quark mass limit. As already explained, 
when |(5i3)_rl| value increase, the CP asymmetry S increase and reach its maximal value at | (^13) jjz, | ~ 0.003 and 
then decrease again. Here we notice two points. First, |((5i3)_rl| is strongly constrained by the inclusive branching 
ratio. In fact, the branching ratio bound is much stronger than the bounds on Ams and sin2/3 m i x . This is because 
b — > dj decay rate is very sensitive to (#13) rl mass insertion due to the enhancement factor nig /mi, as one can check 
it from Eq. (13). On the other hand, there is no such enhancement factor for the B-B mixing process. Second, 
a large CP violation (up to <~ 35%) is still possible even when the only source of CP violating phase is the CKM 
mixing matrix, i.e, when C'^ (nib) is real. The branching ratio Br(B — > X d j) and CP asymmetry S provide the 
complimentary information on (613) rl- 

Figure 3 and 4 are the counterparts of Fig. 1 and 2 respectively, when (5\3)rr dominates. The black region and 
the grey (green) region are defined as in Fig. 1 and 2. In Fig. 3, | (<5i3)j?j?| = 0.03 is assumed while the phase ip is 
varied from to 2ir. We find that CP asymmetry S reach ~ ±5% at best. This small CP asymmetry comes from 
the fact that \C' 7 \ is small, compared to |CV|. The CP asymmetry increases a little as we increase |(#i3).r.r|. But we 
cannot increase |(5i3).r.r| enough to produce large CP asymmetries (say, larger than 10%). This is because the Arris 
and sin2/3 m i x constraints are stronger than the branching ratio bounds in the (613) rr dominating case. (Note that 
the opposite is true for the (Si3)rr dominating case.) The phase <p is set to be zero in Fig. 4, while varying |(5i3)_r.r|. 
We can see that the CP asymmetry is so small that it might be difficult to measure it in the near future. 

In short, large CP asymmetry is possible for the (813) rl dominating case. On the other hand, only small CP 
symmetry is allowed for the (813) rr dominating case. The difference of these two cases is due to the enhancement 
factor m g /nib in the case of (613) rl mass insertion. 

V. CONCLUDING REMARKS 

If we observe a sizable CP asymmetry in B° — > p°j decay, it will be a clear evidence of the new physics beyond 
the SM. Although it is hardly expected that the time dependent CP asymmetry of _B° — > p°j will be measured in 
the present -B-factory, it will be achieved in the next generation of B-factory with about 100 times more B mesons 
produced. Due to the agreement of the SM prediction with the present AniR data and the CP asymmetry in 
B — > J/tpK decay, we favor the new physics which contributes less to the B-B mixing but has a strong effect on the 
b — > d-y penguin diagram. In this paper, we showed that the RL mass insertion of squark mixing of the MSSM can 
produce a large CP asymmetry of B° — > p°-f decay process. 
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FIG. 1: The time-dependent CP asymmetry S as a function of the phase of (5i3)rl- \(8i3)rl \ = 0.001 is assumed. The black 
region denotes allowed points while grey (green) region excluded points by the inclusive b — > d^y branching ratio bound. 




FIG. 2: The time-dependent CP asymmetry S as a function of |(<5i3)iiz,|. The phase of (813) rl is assumed to be 0. The black 
region denotes allowed points while grey (green) region excluded points by the inclusive b — > drf branching ratio bound. 
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FIG. 3: The time-dependent CP asymmetry S as a function of the phase of (8i3)rr. | (<5i3 ) _rj? = 0.03 is assumed. The black 
region denotes allowed points while grey (green) region excluded points by the inclusive b — > drf branching ratio bound. 
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FIG. 4: The time-dependent CP asymmetry S as a function of |(5i3)jjjs|- The phase of (Si3)rr is assumed to be 0. The black 
region denotes allowed points while grey (green) region excluded points by the inclusive b —* dr/ branching ratio bound. 



